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We extend our statistical mechanical theory of the glass transition from examples consisting of 
point particles to molecular liquids with internal degrees of freedom. As before, the fundamental 
assertion is that super-cooled liquids are ergodic, although becoming very viscous at lower temper- 
atures, and are therefore describable in principle by statistical mechanics. The theory is based on 
analyzing the local neighborhoods of each molecule, and a statistical mechanical weight is assigned 
to every possible local organization. This results in an approximate theory that is in very good 
agreement with simulations regarding both thermodynamical and dynamical properties. 



I. INTRODUCTION 

The glass transition is a rather dramatic phenomenon 
in which the viscosity of super-cooled liquids shoots up 
some 14 or 15 orders of magnitude over a small range of 
temperature In fact, the phenomenon is so extreme 
that many authors fitted their measured relaxation time 
to the Vogel-Fulcher formula which predicts a total ar- 
rest of any relaxation at some finite temperature. Recent 
thinking does not support the existence of a finite tem- 
perature singularity. In [2| it was argued that in simple 
glass-forming models the system remains ergodic at all 
temperatures, for any finite number of particles. In Q it 
was shown that even at T = amorphous solids are not 
purely elastic, having some viscous plastic response still 
there even in the thermodynamic limit. In a series of pa- 
pers it was shown that the subtle structural changes that 
occur in simple models of super-cooled liquids are quite 
well described by an approximate statistical mechanics 
that is based on a reasonable up-scaling method that 
defines discrete partition sums involving quasi-species of 
particles and their neighbors All these indicate that 
it may be fruitful to proceed under the assumption that 
the glass transition does not involve the breaking of er- 
godicity, and that therefore statistical mechanics can be 
applied, in fact making use of the extremely slow relax- 
ation to provide a quasi-equilibrium theory in spite of 
the very slow aging [5| . The aim of this paper is to ex- 
tend this approach to the glass transition in molecular 
liquids. These provide additional challenges that do not 
exist in simple models involving point particles, mainly 
due to the existence of of supplementary degrees of free- 
dom and of more complex interactions between the basic 
constituents. We will show, using the example of OTP 
(ortho-terphenyl) which is a very well known glass for- 
mer whose glass transition is extensively studied, that 
the program discussed above continues to make sense, 
providing very needed insights into the phenomenology 
of molecular glass formers. 

The structure of the paper is as follows: In Sect. [ll] 
we introduce the OTP system, explain the model used, 
and present the simulation results that need to be under- 
stood. In Sect, mil the statistical mechanical theory is 
introduced and explained, together with extensive com- 



parisons of the predictions of the theory to the results 
of simulations. In Sect. IIVI we address in particular the 
temperature dependence of the viscosity of super-cooled 
OTP in light of the proposed theory. Finally, Sect. |V] 
is dedicated to a summary and a discussion of the main 
results of this paper. 



II. OTP AND ITS THERMODYNAMICS 

A. The system 

One of the best studied glass-forming system is ortho- 
terphenyl (1,2-diphenylbenzene), denoted for brevity 
OTP. In general, terphenyls (molecular formula Ci8i?i4 
with molecular mass M = 230.31(?/moZ) are a group 
of closely-related aromatic hydrocarbons. They consist 
of a central benzene ring substituted with two phenyl 
groups. There are three isomers, ortho-, meta-, and para- 
terphenyl. 

The structure of the OTP molecule was studied in Refs. 
by different experimental methods and is shown in 
Fig. [H The melting temperature of OTP is T„, = 330K 
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FIG. 1: Color online: The structure of the ortho-terphenyl 
molecule (atomic coordinates from @|). 

and the boiling temperature is Tf, = 605 A' The 
"glass transition temperature" Tg defined by volumet- 
ric and thermal measurements (cf. Figs. I3I5I below) is 
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near 243K [lT| (This temperature depends on the ther- 
mal protocol by which the glass is prepared [l^). 




FIG. 2: Color online: 
terphenyl molecule. 



The planar LW-model of the ortho- 



4>{^cut) = and the first derivative (j)'{rcut) = 0. For the 
cut-off distance r^ut = 1.2616 nm the parameter values 
are a ^ 0.461 kJ/mol, b = -0.313 kJ/(mol nm) 0. 

In the frame of the L-W model of OTP (see Fig. [2]) 
the sites are placed at the vertices of a rigid isosceles 
triangle. The short sides of the triangle are of length a 
and the angle between them is 75°. The position of a site 
a (a = 1, 2, 3) in molecule i is defined by 



p per 



(3) 



where i?^" denotes the position of the center of mass of 
molecule i. In the L-W model OTP molecules are treated 
as rigid bodies; each molecule has six (three rotations and 
three translations) degrees of freedom. Three of them are 
taken to be the Cartesian coordinates of the molecule 
center of mass (define the vector i?J™) and other three 
are the Euler angles specifying the rotational position 
about the center of mass (defined by the vectors ri^a)- 

The energy of interaction between two molecules is 
given by 



B. The model employed in simulations 

A simple coarse grained rigid model of OTP, known as 
the Lewis- Wahnstrom (L-W) model, was introduced for 
the first time in Refs. [l3,E3l- Some intensive studies 
of the properties of this model were achieved by molec- 
ular dynamics (MD) simulations [l3l - [l8| . More realistic 
representation of the OTP molecule consists in taking 
into account internal degrees of freedom and a number 
of models of such kind were developed and used in MD 
simulations [Mill. 

For the sake of efficient numerical simulations the full 
resolution of the structure of the OTP molecule using 
a sufiiciently large number of molecules is beyond our 
computer capabilities. One needs to invoke a reasonable 
model that can retain the essential aspects of the molec- 
ular structure while allowing reasonable computations. 
Therefore we choose the simple Lewis- Wahnstrom model 
and use it for Monte Carlo (MC) simulations. The OTP 
molecule is represented by a three-site complex, each site 
playing the role of the whole benzene ring (see Fig. [Ij . 
The interaction between two sites of different molecules 
is represented by a Lennard- Jones potential 



4>L.j{r) = 4e 



(1) 



where r is the distance between sites. For the sake of 
simulation speed we employ a smooth short range version 
of the interaction potential as proposed in [l6j 



(l>Lj{r) + a + br if r < rcut 
if r > rcut 



(2) 



Here, e = 5.276 kJ/mol and cr = 0.483 nm. The pa- 
rameters a and b are chosen to guarantee the conditions 



^^0(|i?,;,„-4^|). 



(4) 



a = l 13=1 



C. Simulations details 

We simulated iV = 64 and TV = 512 L-W OTP 
molecules in a cubic box with periodical boundary con- 
ditions in the frame of the Monte Carlo method j23 ). 
The trial displacement of a molecule consists of usual 
random translation and random rotation (for details of 
rotations using Euler angles see [Isj). The temperature 
dependence of the particle number density was estimated 
by simulations in an isothermal-isobaric (NPT) ensem- 
ble [IH with the pressure value fixed at P = 1 bar. 
At each temperature the density was used in a canonical 
(NVT) ensemble Monte Carlo simulation. Having results 
in both ensembles enables us below to compare the spe- 
cific heats at constant pressure and at constant volume 
to experiments. After short equilibration the average val- 
ues of the quantities of interest in both ensembles were 
measured during 10^ sweeps in the case of a system with 
iV = 64 and 5 x 10^ for N = 512 molecules. The ac- 
ceptance rate was chosen to be 30%. A final configura- 
tion of MC simulations at constant pressure for = 64 
particles in the simulation cell was used to create an ini- 
tial configuration for simulation with the larger particle 
number A^ = 252 (by using eight small cells to initiate 
one large cell). Final configurations of MC simulations 
at constant pressure were used as initial configurations 
for simulations at constant volume. The following di- 
mensionless units are used: the dimensionless length is 
r* ~ r/a, the particle number density p* = pcr^, where 
p = N/V (V is the system volume), the temperature 
T* ~ ksT/e, where fc^ is the Boltzmann constant, the 
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FIG. 3: Color online: The particle number density as a func- 
tion of the temperature from MC simulations at constant pres- 
sure. Blue triangles and red dots refer to A'^ = 64 and A'^ = 512 
respectively. The solid black line shows the fit to the results of 
the experimental measurements for OTP ([23] and references 
there). Experimental values for the crystal state of OTP [2^ 
are shown for comparison as pink squares. 



energy U* — U/e, the pressure P* — Pa^ /e and the spe- 
cific heat C* = C /{Nks)- In these dimensionless units 



r;, = 0.58, 



1.07 and T* is near 0.43 . 



D. Simulation results. 



The temperature dependence of the particle num- 
ber density at constant pressure is shown in Fig. O 
Obviously, we can not expect that the simple rigid molec- 
ular model of the OTP molecule should reproduce quan- 
titatively the experimental data. Nevertheless, the quali- 
tative behavior of the model is similar to the real system; 
the jump in the slope associated with the glass transition 
[m occurs at the temperature Tg = 0.4, close to the ex- 
perimental value. We have not studied the dependence 
of Tg on the thermal protocol of the sample preparation 
in our simulations (see, e.g., [ll]). The results of Fig. [3] 
are independent of the particle number in the simulation 
cell. This is an indication of negligible finite size effects. 

Fig. m shows simulation results of the potential energy 
change with the temperature at the same conditions as 
Fig. [3] The jump in the slope occurs at the same tem- 
perature Tg = 0.4; again the data are not sensitive to the 
particle number in the simulation cell. 

The specific heat at constant volume, Cy, and constant 
pressure, Cp, are defined by 
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FIG. 4: Color online: The potential energy of OTP model as a 
function of the temperature from MC simulations at constant 
volume. 



where v is the number of degrees of freedom. The 'fluc- 
tuation' part of the specific heat is defined at constant 
volume by 



(6) 



where the potential energy U* is measured in the N-V-T 
ensemble and 



*\2 



NT* 



(7) 



where H* ~ U* —PV is measured in the N-P-T ensemble. 
Simulation results for OTP are displayed in Fig. m 
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FIG. 5: Color online: The temperature dependence of the 
specific heat in the OTP model at constant pressure and con- 
stant volume, such that the volume agrees with the constant 
pressure values (see Fig. [3]). 

Both specific heats show a drop near the temperature 
Tg = 0.4 which is an indicator of the glass transition. 
In experiments one observes qualitative changes in Cp as 
a function of the scanning speed (soj . Our simulations 
are qualitatively similar to what is observed in experi- 
ments with high scanning speed [ll[. With decreasing 
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of the scanning speed the maximum of Cp near the glass 
transition point disappears [l^l ■ We are not aware of ex- 
perimental results on c„ for low temperature OTP; how- 
ever, similar behavior to ours near the glass transition 
point was observed in simulations of binary mixtures (see 

in Hi). 

In order to look for possible structural peculiarities of 
OTP above Tg we studied the radial pair distribution 
function of the molecular center of mass at different tem- 
peratures. Results for three temperatures are shown in 
Fig. [6l At high temperatures this function exhibits an 
ordinary structure typical to liquids with spherical par- 
ticles. With decreasing the temperature the position of 
the first minimum is unchanged, however, the first peak 
displays a development of an internal structure with dips 
at low temperature. Similar behavior was detected in 
Molecular dynamics simulations [l^. This can be at- 
tributed to the consequences of nonsphcrical intcrmolcc- 
ular interactions; at low temperatures and high densities 
the system manifests local orientational ordering. Similar 
evolution of the pair distribution function was observed 
in simulations of hard spherocylinders and it was shown 
that at high packing fraction there exist short range ori- 
entational l33i. 
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FIG. 6: Color online: The radial pair distribution function of 
the OTP center of mass. 

Finally we considered the coordination number defined 
for three dimensional systems as [s^l 

(N) = inp J g{r)r^dr. (8) 



Here r,„i„ is the position of the first minimum of the pair 
distribution function which is one of the possible measure 
of the range short order. The temperature dependence of 
the coordination number is important in thermodynamic 
models of the liquid state [ssl ] . The position of r,„i„ is in- 
dependent of temperature (see Fig. (5]) , therefore, another 
way to estimate the coordination number is to count the 
average number of neighbors inside the sphere of the ra- 
dius r„iin (numerical result coincides with the direct in- 
tegration in Eq. ([5])). The result is displayed in Fig. [T] 
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FIG. 7: Color online: The temperature dependence of the 
coordination number. 



In the liquid state the temperature dependence is typical 
(see UHl)- In the vicinity of Tg and below this value the 
slope is changed. 



III. STATISTICAL MECHANICS OF THE 
GLASS TRANSITION 

A. Statistical geometry approach. 

The structural properties of condensed matter can be 
investigated using scattering experiments (X-rays, elec- 
trons or neutrons scattering, see, e.g., (36|). The angu- 
lar dependence of the scattering intensity is defined by 
the structure factor which is related to the pair distribu- 
tion function. The scattering experiments on liquids and 
amorphous solids exhibit the presence of the short order 
in contrast to periodic crystals with long range order. If 
in a many-body system energy is presented only by two- 
body contributions the estimation of average thermody- 
namics values can be reduced to integrals involving the 
pair distribution function. Unfortunately, usually it is 
impossible to extract the pair distribution function from 
experimental data with the desired accuracy. This func- 
tion is an average property and in the case of a multi- 
component systems even the interpretation of the struc- 
ture of a system is difficult. Looking for example at Fig. 
|6] one sees that the radial pair distribution function is 
not a sufficiently revealing tool to allow a comprehensive 
discussion of the glass transition. As a replacement pos- 
sible tool we turn now to the analysis of the structure 
of the super-cooled liquid in terms of local structures or 
quasi-species as defined below. 

The first suggestion to characterize the state of liq- 
uids by local structures was offered by Bernal in (sj in 
the context of a model of hard balls. The idea is to 
invoke a geometrical approach to define the liquid struc- 
ture. In contrast to a single component crystal with long- 
range order where all particles have the same number of 
neighbors, a liquid is random such that the number of 
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neighbors is not fixed, and can vary quite significantly. 
Defining the concentration of central particles with N 
nearest neighbors by Cat with N G [-/V,nin, ^max], one un- 
derstands that these concentrations will depend on the 
temperature and pressure, defining the 'state' of the liq- 
uid. This approach, which is referred to as 'statistical 
geometry', docs not enjoy an exact method to calculate 
the distribution of concentrations Cn. On the other hand 
this distribution is readily obtained in computer simu- 
lations. In a series of papers (see for example (3] and 
references therein) it was shown that for glass formers 
made of point particles with soft interactions one can 
construct an approximate statistical mechanics that al- 
lows a reasonably accurate description of the evolution 
of the concentrations Ctv as a function of the temper- 
ature, with interesting implications for the study of the 
glass transition. In this paper we extend this approach to 
molecular liquids as explained below. The advantage of 
this approach (see for example [13 ) compared to say, ap- 
proaches based on the pair distribution function, is that 
once the statistical mechanics is set up one can compute 
thermodynamic properties in the usual way that is found 
in any textbook on statistical physics. 



the neighbors with the central molecule. The interaction 
between the neighbors themselves will be estimated be- 
low as a 'crowding effect' that will be characterized by a 
single additional parameter, following ideas prevalent in 
polymer physics. We turn now to a making these ideas 
concrete. 



C. Statistical Mechanics of Quasi-species 
Concentrations 

We begin by assuming that the neighboring molecules 
arc not interacting with each other. In this case every 
neighboring molecule can have an energy taken from m 
possible discrete values ei,-- - Denote by A^i the 

number of ways that we can fit molecules with energy 
£1, consistent with N2 ways to fit molecules with energy 
£2 etc. such that in total Ni + N2 + ■ ■ ■ + N,n = N 
where N < Nmax and Nmax is the maximum number 
of possible nearest neighbors. Then the number of such 
configurations is 



B. Quasi-species in Supercooled Molecular Liquids 

The crucial step in constructing an approximate the- 
ory is the identification of appropriate quasi-species. In 
the case of molecular liquids this task is not immedi- 
ately obvious. As said above, we want to identify quasi- 
species which are made of a central molecule and a num- 
ber of nearest neighbors, and these should have well de- 
fined configurations with well defined energies. The cru- 
cial assumption is that a given quasi-species interacts 
weakly with other quasi-species compared to the inter- 
action of the central molecule with its nearest neigh- 
bors. As attractive as this this idea may sound, there 
are practical problems which appear if the number of 
such quasi-species required to describe the liquid con- 
figurations is not small. For example a non spherically 
symmetric molecule can take up many orientations rel- 
ative to a given central molecule. Each orientation will 
result in a different energy of interaction. Take for ex- 
ample a situation in which the central molecule can have 
up to Nmax neighbors, and each neighbor molecule can 



o(iVi,...,iv„0 



* m.n.T. • 



N I 



As the energy of such a configuration is 

m 

C/(iVi,...,iV„) = ^iVfc£fc, 
fc=l 

the Boltzmann weight of such a configuration is 

m 



(9) 



(10) 



feCfc 



(11) 

and consequently the total weight of any iV nearest neigh- 
bors is 



be in two orientations. Then we already have 2^""''^ pos- W{N) = ^ • • • ^ W{Ni , . . . , Nm)Sj^ j2]^ ^ Nk i (12) 



sible configurations, each in principle having a different 
energy. Thus if Nmax = 35 we might have w 3 x 10^" 
quasi-species-a ludicrous number that destroys the en- 
tire concept of quasi-species. To simplify things towards 
a workable theory (that will have to be validated against 
numerical data) we will treat all nearest neighbor sites as 
being identical, and the molecules occupying these sites 
as having m internal degrees of freedom (for example ori- 
entations relative to the central molecule) with energies 
£1, . . . , £to. The total energy of the quasi-species will be 
estimated as the sum of the energies of interaction of 



where dij is the Kronecker delta. Using the identity 
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Eq. (IT2|) can be written explicitly as 

N ' 

^ ^ max- 
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(14) 



Next wc need to consider the interaction between 
neighboring molecules. Clearly, the assumption that the 
energy of a configuration is defined by Eq. dTU]) is cor- 
rect only for N <C Nmax- The most obvious physical 
effect that destroys this assumption is the soft-core repul- 
sive interaction that will exist between nearest neighbor 
molecule, especially as N approaches Nmax- To include 
such repulsive forces exactly is extremely cumbersome. 
In principle, if all nearest neighbor sites are not iden- 
tical, the repulsive energies can vary between different 
pairs of nearest neighbor molecules. Even if these ener- 
gies are identical the total repulsive energy will depend 
on the exact placement of N nearest neighbor molecules 
in the available Nmax sites. Such considerations go be- 
yond our desire to simplify the model as much as pos- 
sible. We propose therefore a mean-field approximation 
in which the repulsive energy is described by only one 
new parameter. The reasoning is similar to that used in 
the polymer physics [s^, [s^ . Suppose that there exists a 
soft core repulsive energy J between each pair of nearest 
neighbor molecules in the first shell surrounding a cen- 
tral molecule. The probability that a nearest neighbor 
site is occupied is N/Nmax and therefore there will exist 
an additional term to the energy given by Eq. (|10l) due 
to this repulsive contribution. 
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which is quadratic in N. Therefore, Eq. ([H]) is modified 

by the multiplication with the term £xp{ — Y^2W' — ) 
and the concentration of quasi particles with N nearest 
neighbours reads 



C 



N I 



cxp(ln S{T)N 
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where the partition function is 
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* m.n.T. • 



Z{T) 



(16) 
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Z{T) = E 



N=l 



{Nmax - N)\N\ 



X exp \xiS{T)N - 



_J_ N^ 

keT 2Nm 



(17) 



D. Numerical results. 

The most precise way to define 'nearest neighbors' is 
via the construction of a Voronoi tessellation (of polyhe- 
dra or polygons in three- or two- dimensional space cor- 
respondingly) [iOl . In the frame of the Bernal approach 
this method was used in MC simulations to st udy the 
statistical geometry of the Lennard- Jones liquid (4lll4^. 
Unfortunately, the computational cost of a Voronoi con- 
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FIG. 8: Color online: Chekup of Eq. [19] {N„ 
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struction is very high. A simpler definition of 'nearest 
neig hbors' in MC or MD simulations was proposed in 
[43| . The nearest neighbors of a given particle are con- 
sidered as particles with their centers inside the first co- 
ordination shell. It was shown that if the second coor- 
dination shell is excluded the results are not sensitive to 
details of the nearest neighbors definition. Therefore, we 
counted a particle as a nearest neighbor of a central one 
if the distance between their centers of mass is less then 
rmin h^ Eq. [51 The result of this calculation is given in 
terms of the temperature dependence of the concentra- 
tion distribution Cn. 



fill 



{N + 1)Cn^ 



{N„ 



N)C_ 



N 



-] =ln5'(r)- 



J 2N +1 

keT 2Nmax 



(19) 



This suggests that if In [jj^^-^-^^^zw)U^] plotted against 
2N + 1 there should result a straight line of slope 
and intercept InS'(T). Typical simulations 



fcsT 2N,„ 

results and comparison with the prediction of Eq. [19] are 
shown in Fig. \8\ The slope and the intercept calculated 
by a least squares fit for different temperatures are plot- 
ted in Fig. M 

From Eq. p8| the function S{T) is given as a sum 
of exponentials. It turns out that it is very difficult 
to fit S{T) as a sum of a small number of exponen- 
tials with temperature independent energies e^. This 
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FIG. 9: Color online; Upper panel; The function S{T*) de- 
fined as the intercept in lines from Fig. [H] Lower panel; The 
parameter Jdefined as the slope in lines from Fig. [S] . 



may result from the fact that there are many such en- 
ergies, or, as these are in fact phenomenological param- 
eters, they are not at all guaranteed to be temperature 
independent, since they encapsulate information about 
the complicated interactions in the system. The same 
is true for the parameter J. Luckily, in the frame of 
the simple model proposed here we do not need the in- 
dividual parameters since they always appear in the 
same combination in the form of the single function S{T) 
(see Eg. [TB|. This function together with J{T), both of 
which can be estimated from the simulations, are all that 
we need to determine the temperature dependence of the 
concentrations Cm{T) which are the crucial observables 
in the present approach. 

Having determined ^(r) and J(T) from the numerics 
(cf. Fig. in]), the model prediction of the temperature 
dependence of Cm{T) is shown in Fig, 1101 in compar- 
ison with the direct numerical simulation results. The 
agreement appears excellent. It is interesting to com- 
ment here that similar Gaussian looking distributions for 
a given temperature as a function of N was found in all 
the models of glass formers that were studied using a 
similar quasi-species approach Close to the depen- 
dence displayed in Fig. [TU] was observed in simulations 
of water model 
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FIG. 10; Color online: Temperature and dependence of 
Cn{T). Symbols are the direct numerical simulations and 
dashed lines are correspond to the model prediction, Eq. [16] . 



and glass formers once the appropriate quasi-species are 
identified. The difference between one system and the 
other will be encoded in the temperature dependence of 
these Gaussian shapes. To see this we turn in Fig. [11] 
to a comparison of the direct numerical simulations to 
the model predictions for the temperature dependence of 
Cn{T) at chosen values of N. This may be the most 
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44| . We believe that this is a generic fea- 



ture that will be common to a very wide class of liquids 



FIG. 11; Color online; Dependence of the quasi-species con- 
centration on the temperature at fix number of near neighbor 
particles. The dashed lines correspond to the model predic- 
tion Eq. [16] . 

relevant encoding of the subtle changes in the structure 
of the liquid as temperature is changed. We see (quite 
generically again) that some quasi-species decay in their 
concentration when temperature is decreased, some in- 
crease in concentration, and some first increase and then 
decrease to a finite value as T — > 0. The quasi-species 
whose concentration tends to zero are referred to here 
and in previous work as 'liquid-like', since their concen- 
tration is significant only in the high-temperature liquid. 
These are those associated with the highest free energy 
(made from both energy and degeneracy) and therefore 
their concentration declines when the temperature de- 
creases. We will argue that they carry with them many 
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of the signatures of the glass transition. 

Like in previous studies [1^ it is worthwhile to group 
together all the quasi-species into three distinct groups 
according to their qualitative temperature dependence 
(decreasing, increasing, and first increasing and then de- 
creasing). We refer to the resulting model as a 'Doubly 
Coarse Grained' (DCG) description. In the present case 
this means summing the following three groups 



TABLE I: Color online: Partial volumes {v'')l, densities pj, 
and energies {U*)l- 
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35 



N=2 



JV=15 



JV=18 



(20) 

The temperature dependence of these three groups are 
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Ci=Y, Cn{T) ,Cii= J2 Cn{T) , Cm = ^ Cm{T) . 
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FIG. 13: Color online: Volume contributions of the quasi- 
species. 



FIG. 12: Color online: Comparison of MC simulations with 
the model prediction for grouped concentrations (see Eq. I20|) . 

shown in Fig. I12i again comparing the model predic- 
tions to the direct numerical simulation with an obvious 
excellent agreement. One can see that the glass transi- 
tion temperature Tg 0.4 corresponds to the intersec- 
tion point of Cii{T*) and Ciii{T*), at this point C/(T*) 
practically dies out. 

The price of the 'nearest neighbors' definition used by 
us in contrast to the construction of a Voronoi tessella- 
tion is the loss of the information on a volume ascribed 
to a quasi-species. Nevertheless, one can try to use the 
definition of the perfect solution molar volume [46j 

III 



L=I 
III 



(21) 



L=I 



where Wj*^^ = 1/p*. The partial volumes {v*)l were es- 
timated by least squares fits in the whole temperature 
range and are displayed in Table HI The corresponding 
densities = 1/{v*)l are also presented in this table; 
one can see from the comparison with the data of Fig. [3] 
that the value of p}jj is close to the density of the crys- 
tal. Results for the temperature dependence of v1{T*) 
are shown in Fig. 1131 The estimation of the molar vol- 
ume by Eq. (PT|) is compared with MC data. 



In the frame of the approximation of the perfect solu- 
tion the potential energy is given by 



/// 



U* = EC/£(T*) 

L=I 
III 



(22) 



L=I 



The coefficients {U*)l are estimated again by least 
squares fits in the whole temperature range and are dis- 
played in Tab. HI The contributions to the potential en- 
ergy and the comparison of the potential energy from 
simulations and the prediction of Eq. (|22l) are shown in 
Fig. [M]) . The agreement of MC results and predictions 
of Eq. (pij) and Eq. ([^^ means that in spite of the strong 
assumption of additivity in these equations the mixture 
of quasi-species can be effectively considered as a perfect 
solution. 

The enthalpy can be connected with the quasi-species 
concentration using Eq. (|21l) and Eq. (f^^ . Numerical 
differentiation of concentrations given by Eq. ^U\\ yields 
the quasi-species contributions (see Eq. ([5])). Results of 
these calculations and comparison with numerical differ- 
entiation of the enthalpy estimated in MC runs and data 
from Fig. [5] is displayed in Fig. [151 It follows from Fig. 
[T3]-Fig. [15] that the glass transition temperature defined 
by volumetric and thermal measurements can be associ- 
ated with the disappearance of the quasi-species of kind 
I. 
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FIG. 14: Color online: Energy contributions of the quasi- 
species. 




0.2 0.4 0.6 0.8 

T* 



FIG. 15: Color online: Contributions of the quasi-species to 
the heat capacity. 



IV. LOCAL STRUCTURES AND VISCOSITY 
OF OTP 

The theory of the viscosity of gases and liquids saw 
a long and confusing history, especially for super-cooled 
liquids, (irl - lssj . For our purposes we refer to the inter- 
polation formula that was proposed in [ggI - Igs} 

V = Vo.e^^'''- (23) 

Here w/ is the free volume that in general is not easy 
to assign to any physical object. The idea behind this 
interpolation formula is that for particles to move out of 
their immediate cages some free volume needs to open 
up to allow flow. In our context it can be assumed that 
the decline in concentration of Ci{T*) is responsible for 
the increase in viscosity; the free volume in this case is 
defined as v*j = v*j{T*) (see also [H [7|). In order to 
confirm this suggestion it is necessary to consider the 
temperature dependence of the model free volume and 
the measured viscosity from the point of view of Eq. . 

The temperature dependence of the dimensionless vol- 
ume v*j (T) shown in the upper panel of Fig. II 71 above the 



glass transition temperature Tg can be approximated by 




Tg/T 

FIG. 16: Color online: Upper panel: Temperature depen- 
dence of the volume Vj{T). Lower panehTemperature depen- 
dence of the experimental OTP viscosity [7l[ (the value of 7700 
is taken from \7^). 

The linearized temperature dependence of the mea- 
sured OTP viscosity [7l|] is displayed in the lower panel 
of Fig. [ini In the vicinity of the glass transition temper- 
ature it is approximated by 

Inln^ = -2.45 + ^:^, (25) 

where 7700 = 1-25 x 10^'^ poise is taken from [t^I for 
T > 2m.3K. Combination of Eq. ^ and Eq. ([25]) 
yields the free volume expression in the following form 

V = Vo.exp{v:/v*{T)) , (26) 

where v* ~ 2.36. 

In Fig, [17] we present the predictions of the last equa- 
tion (j26p to the experimental measurements of the vis- 
cosity in OTP as a function of T /Tg. The noise in our 
measurement of Ci{T) and vi{T) (cf Figs. 12,13) is re- 
flected in the scatter of the predicted viscosity. If we use 
the model prediction for C/(T) and vi{T) we find the 
solid red line. It follows from this figure that proposed 
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FIG. 17: Color online: Comparison of the experimental mea- 
surements of the OTP viscosity [7l|] with the predictions of 
Eq. (|26|) . Blue rhombi correspond to MC values of concen- 
trations in Eq. H20p . the solid red line represents calculation 
results with concentrations predicted by the model given by 
Eq. (Uni). 

definition of the free volume yields good agreement with 
the experimental data in the whole temperature range 
using the concentration as measured. The model predic- 
tion underestimates the viscosity at low temperatures; 
this must follow from a discrepancy between the model 
prediction and the actual simulational concentrations at 
these temperatures. Nevertheless, the simulation data 
indicates that below the glass transition temperature the 
rate of increase of the viscosity is slowing down. This 
conclusion agrees with the data in [73j . 

V. DISCUSSION. 

The basic assertion of this paper is that super-cooled 
liquids, be them simple or molecular liquids, are ergodic 
systems that arc dcscribablc by statistical mechanics. It 
is quite impossible to write and to solve for the full par- 
tition function that takes into account all the interac- 



tions and the degrees of freedom in these systems. We 
have therefore developed a program in which we offer 
an approximate statistical mechanics that is based on 
considering, for every particle in the system, only the in- 
teractions with the nearest neighbors. This defines our 
'quasi-spccics' which consist of a central particle and its 
N neighbors, where can vary quite a bit. In many 
cases it turns out also advantageous to group the quasi- 
species into groups, such that one group consists of all 
the quasi-species whose concentrations decline when T 
reduces. Usually another group contains all the quasi- 
species whose concentration increase as T decreases, with 
a third group whose concentrations neither increases nor 
decreases We referred to the resulting model as the 
DCG model, in which the concentration of the first group 
acts as an order-parameter that signals the glass transi- 
tion. In other words, since the quasi-species that disap- 
pear when r — >■ are those of highest free energy, we 
consider them as them as liquid-like. In previous pa- 
pers Q we considered the inverse concentration of this 
group as a typical scale measuring the distance between 
the quasi-species raised to the appropriate power. Here 
we opted to compute numerically the volume associated 
with these disappearing objects and considered this as 
the 'free volume' of Eq. (f^. Using this we could es- 
timate the viscosity, and the results were compared to 
experiments in Fig. 1171 It is quite gratifying to see the 
agreement between the predicted viscosity (blue rhombi) 
and the experimental data in the regime below the glass 
transition. We do not have data for T <Tg from the ex- 
periment, but we still have (admittedly quite scattered) 
data for the functions S{T) and J(T) from which we can 
estimate the concentrations of the quasi-species. Using 
these we predict that the viscosity remains finite, even 
though extremely large (note the logarithmic scale in Fig. 
fTTl). even for T — >■ 0. We stress however that what really 
happens at T = cannot be safely deducted from the 
results of this study; it is much better to consider the 
theory of elasticity at T = as can be found in Ref . [3| . 



[1] A. Cavagna, Supercooled liquids for pedestrians. Physics 

Reports 476, 51 (2009). 
[2] J-P. Eckmann and I. Procaccia, Ergodicity and slowing 

down in glass-forming systems with soft potentials: No 

finite-temperature singularities. Phys. Rev. E 78, 011503 

(2008). 

[3] H.G.E. Hentschel, S. Karmakar, E. Lerner and I. Procac- 
cia, Do athermal amorphous solids exist?", Phys. Rev. E 
submitted, also: arXiv:1101.0101 

[4] see for example L. Boue, E. Lerner, I. Procaccia, J. Zyl- 
berg. Predictive statistical mechanics for glass forming 
systems. J. Stat. Mech., PllOlO (2009) and references 
therein. 

[5] D. Dhar and J. L. Lebowitz, Restricted equilibrium en- 
sembles: Exact equation of state of a model glass. Euro. 



Phys. Lett. 92, 20008 (2010. 
[6] C. J. B. Clews, K. Lonsdale, Structure of 1.2- 

Diphenylbenzene (Ci8-ffi4). Proc. R. Soc. Loud. A 

161, 493-504 (1937). 
[7] I. L. Karle, L. O. Brockway, The structure of Biphenyl, 

o-Terphenyl and Tetraphenylene. J. Am. Chem. Soc. 66 

1974-1979 (1944). 
[8] S. Aikawa, Y. Maruyama, Y. Ohashi, Y. Sasada, 1,2- 

Diphenylbenzene (o-Terphenyl). Acta Cryst. B34, 2901- 

2904 (1978). 

[9] G. M. Brown, H. A. Levy, o-Terphenyl by neutron diffrac- 
tion. Acta Cryst. B35, 785-788 (1979). 
[10] W. E. Bachmann, H. T. Clarke, The mechanism of the 
wurtz-fitting reaction. J. Amer. Chem. Soc. 49, 2089- 
2098 (1927). 



11 



R. J. Greet, D. TurnbuU, Glass transition in o-terphenyl. 
J. Chem. Phys. 46, 1243 (1967). [31 
S. S. Chang, A. B. Bestul, Heat capacity and thermody- 
namic properties of o-terphenyl crystal, glass, and liquid. 
J. Chem. Phys. 56, 503-516 (1972). [32 
L. J. Lewis, G. Wahnstrom, Relaxation of a molecular 
glass at intermediate times. Solid State Comm. 86, 295- 
299 (1993). [33 
L. J. Lewis, G. Wahnstrom, Molecular-dynamics study of 
supercooled ori/io-terphenil. Phys. Rev. B 50, 3865-3877 [34 
(1994). 

A. Rinaldi, F. Sciortino, P. Tartagha, Phys. Rev. E 63, 
061210 (2001). [35 
S. Mossa, E. La Nave, H. E. Stanley, C. Donati, F. 
Sciortino, P. Tartaglia, Dynamics and configurational en- 
tropy in the Lewis- Wahnstrom model for supercooled or- [36 
thoterphenyl. Phys. Rev. E 65, 041205 (2002). 
S. H. Chong, F. Sciortino, Structural relaxation in super- 
cooled orthoterphenyl. Phys. Rev. E 69, 051202 (2004). [37 
T. G. Lombardo, P. G. Debenedetti, F. H. Stillinger, 
Computational probes of molecular motion in the Lewis- [38 
Wahnstrom model for ortho-terphenyl. J. Chem. Phys. 
125, 174507 (2006). 

S. R. Kudchadkar, J. M. West, Molecular dynamics sim- [39 
ulations of the glass former ortho-terphenyl. J. Chem. 
Phys. 103, 8566-8576 (1995). [40 
S. Mossa, R. Di Leonardo, G. Ruocco, M. Sampoli, 
Molecular dynamics simulation of the fragile glass-former 
orthoterphenyl: A flexible molecule model. Phys. Rev. E [41 
62, 612-639 (2000). 

S. Mossa, R. Di Leonardo, G. Ruocco, M. Sampoli, 
Molecular dynamics simulation of the fragile glass-former [42 
orthoterphenyl; A flexible molecule model. II. Collective 
dynamics. Phys. Rev. E 64, 021511 (2001). 
J. Ghosh, R. Faller, A comparative molecular simulation [43 
study of the glass former ortho-terphenyl in bulk and 
freestanding films. J. Chem. Phys. 125, 044506 (2006). 
R. J. Berry, D. Rigby, D. Duan, M. Schwartz, Molecular [44 
dynamics study of translation and rotational diffusion 
in liquid ortho-terphenyl. J. Phys. Chem. AllO, 13-19 
(2006). [45 
N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. 
H. Teller, E. Teller, Equation of state calculations by fast 
computing machines. J. Chem. Phys. 21, 1087 (1953). 
J. C. Owicki, H. A. Scheraga, Monte Carlo calculations [46 
in isothermal-isobaric ensemble. 1. liquid water. J. Amer- 
ican Chem. Soc. 99, 7403-7412 (1977). 

W. W. Wood, Monte Carlo Studies of Simple Liq- [47 
uid Models. In: H.N.V. Temperley, J.S. Rowlinson and 
G.S. Rushbrooke, Editors, The physics of simple liquids, 
North-HoUand, Amsterdam, pp. 115-230 (1968). [ 
G. Monaco, D. Fioretto, L. Gomez, G. Ruocco, Glass 
transition and density fluctuations in the fragile glass for- 
mer orthoterphenyl. Phys. Rev. E 63, 061502 (2001) [49 
J. Opdycke, J. P. Dawson, R. K. Clark, M. Dutton, J. 
J. Ewing, H. H. Schmidt, Statistical thermodynamics of 
the polyphenyls. I. Molar volumes and compressibilities [50 
of biphenyl and m-, o-, and p-terphenyl. J. Phys. Chem. 
68, 2385 (1954). [51 
D. N. Pererra, P. Harrowell, Stability and structure of 
supercooled liquid mixture in two dimensions. Phys. Rev. 
E 59, 5721-5743 (1999). [52 
A. ToUe, Neutron scattering studies of the model glass 
former ortho-terphenyl. Rep. Prog. Phys. 64, 1473-1532 [53 



(2001). 

H. G. E. Hentschel, V. Ilyin, I. Procaccia, Nonuniversal- 
ity of the specific heat in glass forming systems. Phys. 
Rev. Lett. 101, 265701 (2008). 

H. G. E. Hentschel, V. Ilyin, I. Procaccia, N. Schupper, 
Theory of specific heat in glass-forming systems. Phys. 
Rev.E 78, 061504 (2008). 

H. Lowen, Brownian dynamics of hard spherocylinders. 
Phys. Rev. E 50,1232-1242 (1994). 

W. Brostow, Radial distribution function peaks and co- 
ordination numbers in liquids and amorphouse solids. 
Chem. Phys. Lett. 49, 285-288 (1977). 
S. Skjold-J0rgensen, P. Rasmussen, A. Fredenslung, On 
the temperature dependence of the UNIQUAC/UNIFAC 
models. Chem. Eng, Sci. 35, 2389-2403 (1980). 
J. M. Ziman, Models of disorder. The theoretical physics 
of homogeneously disordered systems. Cambridge Uni- 
versity Press (1979). 

J. D. Bernal, The structure of liquids. Proc. R. Soc. Lond. 
A280, 299-322 (1964). 

R. Fiirth, A new approach to the statistical thermody- 
namics of liquids. Proc. Roy. Soc. Edinb. A66, 232-251 
(1963). 

P. de Gennes, Scaling concepts in polymer physics. Cor- 
nell University Press (1979) 

G. F. Voronoi, Nouvelles applications des parametres 
continus a la theorie de formes quadratiques. J. reine 
angew. Math. 134, 198-287 (1908). 

J. L. Finney, Random packing and the structure of simple 
liquids. I. The geometry of random close packing. Proc. 
R. Soc. Lond. A 319, 479-493 (1970). 
J. L. Finney, Random packing and the structure of sim- 
ple liquids. II. The molecular geometry of simple liquids. 
Proc. R. Soc. Lond. A 319, 495-507 (1970). 
P. J. Steinhardt, D. R. Nelson, M. Ronchetti, Bond- 
orientational order in liquids and glasses. Phys. Rev. B 
28, 784-805 (1983). 

S. Merchant, J. K. Shah, D. Asthagiri, Water coordi- 
nation structures and the excess free energy of liquids. 
.arXiv:1101.1076y l [physics.chem-ph] (2011). 
V. Ilyin, E. Lerner, Ting-Shek Lo, I. Procaccia, Statis- 
tical mechanics of the glass transition in one-component 
liquids with an anisotropic potential. Phys. Rev. Lett. 
99, 135702 (2007). 

I. Prigogin, The molecular theory of solutions, avec A. 
Bellemans et V. Mathot; North-Holland Publ. Company, 
Amsterdam, 1957. 

J. O. Hirschfelder, C. F. Curtiss, R. B. Bird, Molecu- 
lar Theory of Gases and Liquids. John Wiley, New York 
(1954). 

A. J. Batschinski, Untersuchungen fiber die innere Rei- 
bung der Fluessigkeiten, Z. Physik Chem. 84, 643-706 
(1913) 

D. B. Macleod , On a relation between the viscosity of a 
liquid and its coefficient of expansion. Trans. Farad. Soc. 
19, 6-16 (1923). 

J. H. Hildebrand, Motions of molecules in liquids: viscos- 
ity and diffusivity. Science 174, 490-493 (1971). 
L. D. Eicher, B. J. Zwolinski, Limitations of the 
Hildebrand-Batschinski shear viscosity equation. Science 
177, 369 (1972). 

E. N. da C. Andrade, The Viscosity of Liquids, Nature 
125, 309-310 (1930). 

H. Vogel, Das Temperaturabhaungigkeitsgesetz der 



12 



Viskositat von Fliissigkeiten. Phys. Zeit. 22, 645-646 
(1921). 

[54] G. S. Fulcher, Analysis of recent measurements of the 
viscosity of glasses. J. Am. Ceram. Soc. 8, 339-355 (1925). 

[55] G. Tamman, Glasses as supercooled liquids. J. Soc. Glass 
Technol. 9, 166-185 (1925). 

[56] G. Tamman, W. Hesse, Die Abhangigkeit der Viscositat 
von der Temperatur bei unterkiihlten Fliisigkeiten. Z. 
anorg. u. allgem. Chem. 156, 245-257 (1926). 

[57] H. Eyring, Viscosity, plasticity, and diffusion as exam- 
ples of absolute reaction rates. J. Chem. Phys. 4, 283-291 
(1936). 

[58] F. Mallamace, C. Branca, C. Corsaro, N. Leone, J. 
Spooren, S-H. Chen, H. E. Stanley, Transport properties 
of glass-forming liquids suggest that dynamic crossover 
temperature is an important as the glass transition tem- 
perature. PNAS 107, 22457-22462(2010). 

[59] G. Adam, J. H. Gibbs, On the temperature dependence 
of cooperative relaxation properties in glass-forming liq- 
uids. J. Chem Phys. 43, 139- (1965). 

[60] R. Richert, C. A. Angell, Dinamics of glass-forming liq- 
uids. Y. On the link between molecular dynamics and 
configurational entropy. J. Chem. Phys. 108, 9016-9026 
(1998). 

[61] A. V. Granato, A derivation of the Vogel-Fulcher- 
Tamman relation for supercooled liquids. J. Non-cryst. 
Solids 357, 334-338 (2011). 

[62] I. Avramov, Viscosity in disordered media. J. Non-cryst. 
Solids 351, 3163-3173 (2005). 

[63] T. Hecksher, A. I. Nielsen, N. B. Olsen, J. C. Dyre, Little 
evidence for dynamic divergences in ultraviscous molec- 
ular liquidss. Nature Phys. 4, 737-741 (2008). 

[64] O. N. Senkov, D. B. Miracle, Descriprtion of the fragile 
behavior of glass-forming liquids with the use of experi- 



mentally accessible parameters. J. Non-cryst. Solids 355, 
2596-2603 (2009). 
[65] I. Avramov, Interrelation between the parameters of 
equations of viscous flow and chemical composition of 
glassforming melts. J. Non-cryst. Solids 357, 391-396 
(2011). 

[66] A. K. Doolittle, Studies in Newtonian flow. II. The de- 
pendence of the viscosity of liquids on free-space. J. Appl. 
Phys. 22, 1471-1475 (1951). 

[67] A. K. Doolittle, D. B. Doolittle, Studies in Newtonian 
flow. Y. Further verification of the free-space viscosity 
equation. J. Appl. Phys. 28, 901-905 (1957). 

[68] M. H. Cohen, D. TurnbuU, Molecular transport in liquids 
and glasses. J. Chem. Phys. 31, 1164-1169 (1959). 

[69] E. Aharonov, E. Bouchbinder, H. G. E. Hentschel, V. 
Ilyin, N. Makedonska, I. Procaccia, N. Schupper, Di- 
rect identification of the glass transition: Growing length 
scale and the onset of plasticity. EPL 77, 56002 (2007). 

[70] H. G. E. Hentschel, V. Ilyin, N. Makedonska, I. Procac- 
cia, N. Schupper, Statistical mechanics of the glass tran- 
sition as revealed by a Voronoi tesselation. Phys. Rev. E 
75, 050404(R) (2007). 

[71] D. J. Plazek, C. A. Beroa, I.-C. Chay, The recoverable 
compliance of amorphous materials. J. non-cryst. solids 
172-174, 181-190 (1994). 

[72] M.T. Cicerone, M.D. Ediger, Photobleaching technique 
for measuring ultra-slow reorientation near and below 
the glass transition; the tetracene/o-terphenyl system. 
J. Phys. Chem. 97, 10489-10497 (1993). 

[73] H. Kobayashi, H. Takahashi, Y. Hiki, Viscosity of glasses 
near and below the glass transition temperature. J. Appl. 
Phys. 88, 3776-3778 (2000). 



